o 



(N 



Shearlet frames with short support 

Song Li*and Yi Shen''^ 



Abstract 



Compactly supported shearlets have been studied in both theory and apphcations. In this 
paper, we construct symmetric compactly supported shearlet systems based on pseudo splines 
of type II. Specially, using B-splines, we construct shearlet frame having explicit analytical 
forms which is important for applications. The shearlet systems based on B-splines also provide 
' optimally sparse approximation within cartoon-liked image. 
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1 Introduction 

> , 

. Cartoon-liked image are 2-dimensional functions that are except for discontinuities along 

curves tlj. To find optimally sparse representations of cartoon-like image, several variations of 
Tj- ■ the wavelet scheme have been proposed, such as curvelets [1] and contourlets [8]. Shearlets frame 

developed in [T7j is the first multiscale directional system which also provides almost optimally 
sparse approximation with cartoon-like images. However, these studies are only concerned band 
limited generator. Very recently, Kutyniok and Lim presented a complete proof of (almost) opti- 
mally sparse approximations of cartoon liked images by using shearlet systems which are generated 
by compactly supported shearlets under some weak moment conditions [16] . They also constructed 
^ ' a class of compactly supported shearlet frames based on pseudo splines of type I [15]. Hence, 

■ excellent spatial localization is achieved. But the shearlet frame still have two disadvantages: 

• The shearlet is not symmetric or anti-symmetric; 

• The shearlets do not have explicit analytical forms in spatial domain. 

These drawbacks motivate us to consider constructing shearlet frame using B-splines. B-splines 
had a significant impact on the development of the theory of the wavelet analysis. They yield the 
only wavelets that have explicit analytical forms. All other wavelet bases are defined indirectly 
through an infinite product in Fourier domain [5l |6] . As Daubechies pointed out in [6] , except the 
Haar wavelet function, there is no compactly supported real-valued symmetric orthonormal wavelet 
basis in L2(M). However, it is much easier and more flexible to construct and design compactly 



'Department of Mathematics, Zhejiang University Hangzhou, 310027, China 

^Corresponding author: syll33@163.com, Department of Mathematics, Zhejiang University Hangzhou, 310027, 
China 



1 



supported wavelet frames or Riesz bases than orthonormal wavelet bases. For example, from any 
B-spline function of order m, one can construct a symmetric tight wavelet frame with m generators 
[18j . Tight wavelet frame from B-splines with high vanishing moments were considered in j2i |7J . 
The compactly supported Riesz wavelets generated from B-splines were first constructed in [3]. 
The shortest supported Riesz wavelet with m vanishing moments from B-spline of order m were 
constructed in [12]. The compactly supported wavelet bases from B-splines for Sobolev spaces were 
investigated in [HI [13] . 

The rest of the paper is organized as follows. In Section [2l we construct shearlets based on 
B-splines, then we present some results on the optimally sparse approximations of cartoon- like 
images. In Section 3, we investigate the lower bounds and the upper bounds of the pseudo splines 
in Fourier domain. These results not only have their own interests but also have closed relations 
to the shearlets frame bounds. In sections 4, we presents some examples to illustrate our results. 



2 Shearelts based on B-splines 

In this section, we first introduce the definitions of shearlet frame and B-splines function. Then we 
construct compactly supported shearlet frame with the generator from B-splines. Finally, we show 
that these compactly supported shearlet systems provide (almost) optimally sparse approximations 
of cartoon-liked images. 

Shearlets are scaled according to a parabolic scaling law encoded in the matrix A2j or A2j, j € Z, 
and exhibit directionality by parameterizing slope encoded in the matrices Sk, k G Z, defined by 

A - f \ . _f 2^/2 \ ^ _f I k\ 

- ^ 2^/^ j ' - ^ ), '^fc - ^ 1 )■ 

Now we define discrete shearlet systems in 2D. Let c = (ci,C2) G (M+)2. For (pjipjip S L'^(M?), the 
cone-adapted 2D discrete shearlet system SH{cj),ip,ip; c) is defined by 

SH{(P, i;, Vi; c) = a) U ^{i^; c) U ^C^i; c), 

where 

ci) = {(/)(■- cm) : m E Z^}, 
^{iP;c) = {2-4^ijiSkA2, ■ -cm) : j > 0, < \2^/^l m G Z^}, 
*(^;c) = {2^^4j{SIA2j ■ -cm) : j>0, \k\ < \2^/'^], m G Z^}. 

We partite the frequency plane into Ci{a) — Ci{a) where 

6 >«, 16/61 < 1} : 1 = 1 

6 > a, 16/61 < 1} : 1 = 2 

6 > -a, 16/61 < 1} : / = 3 

6>-«, 16/61 <1}: 1 = 4. 

and a centered rectangle 

^(«) = {(6,6)K2:||(6,6)l|oo<a}. 



C,(a) 



{(6,6) 
{(6,6) 
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The region Ci UC3 is covered by the frequency support of shearelets in ^{^; c). The region C2 UC4 is 
covered by the frequency support of shearelets in ^{tl^; c). The region TZ is covered by the frequency 
support of (j). Recall {fTjjjg/ form a frame in L^(M^) if there exist constants A,B>0 such that 

^ll/lli<EK/'^*)l'^^ll/ll2 v/gl2(r2). 

The numbers A, B are called frame bounds. If ((/>, ■0, '0; c) is a frame for L^(]R^), we refer to ij: 
and as shearlets. 

We say that (p is a refinable function with mask a(^) if 0(2^) = ^ ^ The 

Fourier transform of a function / € -Li(R'^) is defined as /(^) = f{x)e~'^^'^dx and can be 
naturally extended to tempered distributions. As an important family of refinable functions, B- 
spline functions are useful in applications. B-spline with order m and its mask is defined by 

^(0 = e"^(^^^^)" and a(0 = cos-(C/2), ^ G M, 

where j = when m is even, and j = 1 when m is odd. The B-spline function Em G C™~^(M) is a 
function of piecewise polynomials of degree less than m, vanishes outside the interval [0, m] and is 
symmetric about the point x = m/2. Now we state our first contribution on the shearlet frame in 
this paper. 

Theorem 2.1. LetBj\[ be the B-spline function of order N with the mask a^iO = (2~^)(l+e~*^)^. 
Let Ni,N2 en be such that Ni > N2 > 3. Define a shearlet G L'^{^^) by 

4,(20 = 2-^^e-^Hl - e-«^)^^^,(6)^.fe), C = (6,^2). 

For given < a < it/2, there exits a sampling constant c > such that the shearlet system ^{^; c) 
forms a frame for {/ G L2(M^) : supp/ G Ci{a) U Cs{a)} for C2 < ci < c. 

Proof. The proof is a straightforward consequence of Theorem 13.101 □ 

The cartoon-liked model was first introduced in [1]. The basic idea is to choose a closed boundary 
curve and then fill the interior and exterior parts with functions. For v > 0, the set STAB?'{v) 
is defined to be the set of all i? C [0, 1]^ such that i? is a translate of a set 

{x G : jx| < p{9), X = {\x\,9) in polar coordinates} 

which satisfies \p"{0) < y\, p < po < 1- Then (^^(z^) denotes the set of functions / G L{M?) of the 
form 

/ = /o + flXB, 

where /o, /i G Cq([0, 1]^) and B G STAE?{v). The bandlimited curvelets, contourlets, and shear- 
lets exhibit (almost) optimally sparse approximation with this model. The first complete proof of 
(almost) optimally sparse approximations of cartoon-liked images by using compactly supported 
shearlet frame was given in [16]. Let us now be more precise, and introduce these results. Let 
c > 0, and let (f), ijj, ^ L2(M^) be compactly supported. Suppose that for all ^ = (^1,^2) € 1^^, 
the shearlet ip satisfies 

1^^(61 < ■ min{l, 161"} • min{l, ■ min{l, j^m, (2-1) 
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l4.>(OI<l'.(6)l(i + |^ 



(2.2) 



where a > 5, 7 > 4, /i E Li(M), and Ci is a constant, and suppose that the shearlet ip satisfies 
(j2.ip and (j2.2p with the roles of ^1 and ^2 reversed. Further, suppose that SH(c; 4>, ip) forms a 
frame for L2(M^). Denote {(Ji)i^i = SH{c; (p,tp,ip). Let (cr)jg/ be a dual frame of {cri)i^j. We can 
take the nonlinear A^-terms approximation 



N 



{f,cri)^i 



where ((/, Ci))je/jv largest coefficients {f,(Ji) in magnitude. Then, for any u > 0, the 

shearlet frame SH{c; cf), ijj, tp) provides (almost) optimal sparse approximation of function / € (^'^{v) 
in the sense that there exists some C > such that 

\\f-fN\\l<C-i\ogNf-N-^ 

The condition (j2.1|) and (j2.2p can be viewed as a generalization of a second order directional 
vanishing moment condition, which is crucial for having fast decay of the shearlet coefficients. 
Following the line of ^5j, we obtain the following results which is our second contribution in this 
paper. 

Theorem 2.2. LetBjyj he the B- spline function of order N with the mask a^{^) = {2~^){l+e~^^)^ . 
Let Ni,N2 GN be such that iVi > 5, A^2 > 4 and Ni > N2. Define a shearlet V'l G -L^(M^) by 

M20 = 2-^ie-«Hl - e-«^)^i 5^,(6)^,(^2) ^ = (6,6)- 

Let 4){x) = B]\j^{xi)Bn^{x2) and 'i/'2(a;i> 2:2) = "01(^1) ^2)- Then there exit sampling constant c > 
such that the shearlet system SH{(j),ipi,ip2'-,c) provides (almost) optimally sparse approximations of 
function f € S''^{v) in the sense that there exists some constant C > such that 

\\f-fN\\l<C-{\ogNf-N-\ 

where fjy is the nonlinear N-term approximation obtained by choosing the N largest shearlet coef- 
ficients of f . 

Proof. By Theorem 12.11 there exits a sampling constant c > such that the shearlet system ^'(-0; c) 
forms a frame for {/ G L2(M^) : supp/ G Ci(a) U Cz{a)} with c = (ci, C2) G (1^+)^ and C2 < ci < c. 
With the same argument as in [15], we can prove that SH{c; (j), -0, ip) forms a frame for L2(M) with 
c = (01,02) G (M+)^ and 02 < oi < c. 
Rewrite V'i(C) iii the following 



MO 



i37V2(ei)-(l-e" 



^Af2 (6 



We obtain 



dMO 



96 



2-^1 e- 



< 



(1 



4 



Set = 2-^ie-^«i(l - e-*«i)^i-^2S^2(6). We have 

IM6)l<c"^,(ei)<c"|6 



-N2 



Therefore, € Lii^ 

Let us consider \Bm^ (^2)|- When iV2 is even, Bm^{^2) = (^^^^75^) • Then 



T^'(,^ r sinfe/2) ^^^-^ r cosfe/2) sin(6/2) 



< C\i2 



When iV2 is odd, S^^fe) = e'^^^ (^^Hg^M^^'. Then 



|Siv/(6)l = 



_. /sin(6/2)\^^ ^ sin(6/2) ^^^-^r cos(6/2) sin(g2/2) 

V 6/2 n 6/2 [ 4^2 2^1 



< 



sin 



< C'I6 



^(6/2) 
6/2 

\-N2 



N2 



+ 



No 



- sin(6/2) \^^-^ r cos(g2/2) _ sin(g2/2 



6/2 

constant C such that 



46 



2ei 



We can conclude that there exists a 

< C7min{l, leiT^} • min{l, \^,\-^^} ■ min{l, l^r^^^, 



and 



^6 V i?i . . 

where A^i > 5, A''2 > 4, /i G Li(M). Hence, the shearlet ipi satisfies the conditions (j2.ip and (j2.2p . 
and the shearlet '02 likewise. The theorem is proved n 



/ I/: I \ -A^2 



3 Shearlets based on pseudo splines 

In this section, we first briefly recall a family of refinable function: pseudo splines. After establishing 
some useful lemmas, we investigate the lower bound and the decay of the Fourier transform of 
pseudo splines of type II. Finally, we construct shearlet frame based on pseudo splines. 

Pseudo-splines are compactly supported refinable functions in L2(M). For positive integers 
A^, / € N with / < A^, denote P]\i^i{x) := Yl^j=Q ■ The mask of a pseudo spline of type II 

with order (A^, I) is defined by 

2-^(0 := cos2^(e/2)Piv,Ksm'(e/2)). (3.1) 

The mask of a pseudo spline of type I is defined by 2aN\i{0 •= Hence, fajv^(^) with 

real coefficients is the square root of liajyTKOP using the Riesz lemma. In general, pseudo splines 
of type I are neither symmetric nor antisymmetric. To achieve symmetry, compactly supported 
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complex valued pseudo splines were introduced in [19]. The corresponding pseudo splines can be 
defined in terms of their Fourier transform, i.e. 

oo 

'^{0 = ll^i{'^~'0, A: = 1,2. 

i=i 

An important fact is that 2fljv\z(0 defined by the summation of the first / + 1 terms of the binomial 
expansion of (cos2(^/2) + sin2(^/2))^+' = 1 [71[9] i.e. 

E " ' ^ ') = E ^ ''"''(^/'^ cos^('-^)(e/2). (3.2) 

The pseudo splines with order (A^, 0) are B-splines. For the case / = A^— 1, the pseudo splines of type 
I are orthogonal refinable functions given in [6], and the pseudo splines of type II are interpolatory 
refinable function given in |10] . The other pseudo splines fill in the gap between the B-spline and 
orthogonal refinable functions for type I and B-spline and interpolatory refinable function for type 
II. For positive integers N,l £N with / < A^, let a(^) be the mask of the pseudo splines of type II 
with order {N,l). Then oaTKO can be factorized as 

=cos2^(e/2)|/:(01, ee[-^,vr]. (3.3) 

This shows that pseudo splines is the convolution of a B-spline of some order with a distribution. 
Since £(^) is bounded, £(^) is actually the mask of a refinable distribution. The regularity of cp 
comes from the cos^^(^/2) factor. The distribution part provides some desirable properties for 
(j), such as orthogonality of its shifts. In [9], Dong and Shen gave a regularity analysis of pseudo 
splines of both types. The key to regularity analysis is 

1^1 < f , (3.4) 
|/:(0£(2e)| < |£(y)P, f <|e|<vr. 

To investigate the bounds and the decay of the pseudo splines in Fourier domain, we establish the 
following lemmas which have their own interesting. For simplify, we denote the mask of pseudo 
splines of type II by a(^). 



1^(01 < 



3.1 Lemmas 

Lemma 3.1. For positive integers A", Z G N with I < N, let a(^) be the mask of the pseudo splines 
of type II with order (N, I) and let C{(,) be defined as in i3.3\) . Then 



|a(OI>l-C7i|e| 



21+2 



where Ci = ^%t'+\ ' ' 
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Proof. By (j3.2p . we obtain 

N+l 



j=i+i V / 

^—1^1 \ J / 



j=l+i 
sr^N+l (N+l\ 
l^j=l+l \ j J \c\2l+2 



□ 



Lemma 3.2. For positive integers A^, Z € N with I < N , let a(^) be the mask of the pseudo splines 
of type II with order (N, I). Then 

\m\<^ + c2\^\^ 

/jv-l+j\ 

where C2 = '^^^^ ' ' ■ 
Proof. We obtain 

I 



j=o V J / 



sm 

□ 

Lemma 3.3. For positive integers A^, Z € N with I < N , let a{(,) be the mask of the pseudo splines 
of type II with order (N, I). Let b{£,) = e~*^a(^ + ir), then 

\b{a)\xi-p,-aUaMO < \m < min{l, ' 1^^}, < a < /? < vr. (3.5) 

Proof. \b{(,)\ is a TT shift of the 27r-periodic triangle polynomial a(^). Note that |a(.^)| is increasing 
on [— vr, 0] and decreasing on [0, vr]. It is easy to see that is decreasing on [— vr, 0] and increasing 
on [0,7r]. Hence the left side of ([33]) holds. By ([33]), We have 

m\ = sin^^(e/2)5:(^ + ^)cos2^(e/2)sin2('-^)(e/2) 
^ ^ l?l • 

□ 
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3.2 Regularity 

Now we give the lower bound for the pseudo spHnes of type II in the Fourier domain. 
Theorem 3.4. Let <p be the pseudo-splines of Type II with order {N,l). Let K C [— tt, vr], then 

\m\>CA-xK{i). 

where C4 = Hfeii exp(-Ci2-^o+i|^o 

Proof. Let ^0 = argmax^g^ Since a(^) is decreasing on [0,7r], we have 

for fc > 1 and ^ G We choose sufficiently large ko so that 2~'^Ci|^p+^ < ^ ii ^ £ K and k > ko- 
Using 1 - X > e~2=^ for < x < 1/2, we obtain 

\aiO\ > 1 - Ci2-'^|^|2'+2 > exp(-2'=+iCi|4|2'+2), k > kg. 

Therefore, for ^ G 

fco 00 

fc=l k=ko+l 

ka 00 

> Yl\ai2-%)\ n exp(-2Ci2-'=|e|2^+2) 

k=l k=ko+l 
ko 

> n |a(2-'=eo)| exp(-Ci2-*^o+i|^oP'+') = C4. 
fe=i 

□ 

In the following, we mainly investigate the decay of the pseudo splines of type II in the Fourier 
domain. 



Lemma 3.5. For positive integers N,l (^N with I < N , let a(^) be the mask of the pseudo splines 

1% 



of type II with order {N, I). Let C{^) be the distribution part. Assume that (j^ciO — YYjLi '^(2 "'C); 



then ^ 

MO < exp(C2/3)gig^i|^|i°S2(«i (3.6) 

where q, = sup|^|<, = E;=o C^^^') '^^id q2 = \C{27r/3)\. 

Proof Since < 1 + Czj^p, we have for |^| < 1, 

00 / 00 

l[\C{2-^0\ = expK;iog(|/:(2-^-^|) 

(00 
C2Y,^2-^'\^\ 

< exp(C2/3). 
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Now for |,^[ > 1, we have 



n ^2-^0 = n /(2"^'o 



fe=i 



k=0 



where /(^) = n/=i'^(2 "'O- For the given j^| > 1, there exists a positive integer /cq such that 
2^0 J < < 2(fco+i)J. Denote r/ = 2-(''o+i)-^^, we obtain |r?| < 1. Define qi = sup|g|<^ \C{C)\ and 
52 = |/:(27r/3)|. By ([331), we obtain 



< 9192 



Therefore, 



A;=fco + l 



Moreover, 



j=0 

oo 



i=0 



J=0 



ni^(2~^'^)| <exp(C2/3). 

j=0 



n 1/(2- 

A;=0 



-kJ\i ^ / J-1 



fco+1 



< 9192 9l 92 



fco(J-l) 



S 91 92 (9i 92 



< 9i92"^|C|^°S2(^i^*''~''^'^ 



We conclude that 



□ 



Theorem 3.6. For positive integers N,l G 'N with I < N , let (j) be the pseudo splines of type II 
with order (N, I). Then for any given integer J > 1. 



\m\ < min{l,C3|Cr2^+^°S2(^'''''"'«^)} 
where C3 = 4^ exp(C2/3)gig^~^ 
Proof. Since 



(3.7) 



we obtain |0(^)| < 1. 
It is weU known that 



n°° sin(g/2) 
^cos(2 ^0 = 
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Table 1: Decay rate (3^,1 = 2A^ — k of pseudo splines of type II for 2 < < 9 



(iv,o 


iV = 2 


iV = 3 


iV = 4 


= 5 


= 6 


N = 7 


iV = 8 


iV = 9 


/ = 


4.00000 


6.00000 


8.00000 


10.0000 


12.0000 


14.0000 


16.0000 


18.0000 


I = 1 


2.67807 


4.29956 


6.00000 


7.75207 


9.54057 


11.3561 


13.1927 


15.0458 


I = 2 




3.27208 


4.73321 


6.27890 


7.88626 


9.54057 


11.2318 


12.9530 


I = 3 






3.82507 


5.19506 


6.64465 


8.15608 


9.71691 


11.3181 


I = 4 








4.35316 


5.66363 


7.04717 


8.48992 


9.98156 


I = 5 










4.86449 


6.13261 


7.46770 


8.85865 


I = 6 












5.36349 


6.59988 


7.89780 


I = 7 














5.85310 


7.06473 


I = 8 
















6.33529 



Hence 



\m\ 



By (j3.6p . we have 



2N 



exp(C2/3)gi (72^-1 lel 



, , 1/(.J-1) N 

|log2(9i 92) 



where C3 = 4^ exp(C2/3)gig^ 



Corollary 3.7. For positive integers N,l G'N with I < N, let 
with order {N, I). Then 

m)\<Cmin{l,\^r'''+^}, 

where k = log(P(3/4))/log 2. 

Proof. The proof is a straightforward consequence of Theorem 



□ 

be the pseudo splines of Type II 

(3.8) 

□ 



Remark 3.8. The above results were first proved by Dong and Shen in f^. The decay of the 
Fourier transform is optimal. We give estimation 1^3. 7\ ) with explicit constant which is important 
for the frame upbound of the shearlet frame. Table 1 gives the decay rate of the Fourier transform 
of pseudo splines of Type II with order {N, I) for 2 < N < 9 and < I < N . The decay rate of the 
Fourier transform of pseudo splines of Type II with order {N,l) is /37v,«/2. 



3.3 Shearlet Frames 

In the last subsection, we construct the shearlet frames based on pseudo splines. We give a weaker 
condition for constructing shearlet frames than Theorem 4.9 in ^15j . 
For function V', V' € L'^{R'^), we define G : x ^ R by 

e(C,u;) = |<^(e)||0(e + ^)| + ei(e,a;) + 62(^,0;), 
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where 
and 

Q2{C,u) = Y, \^{SkA2-.mkSkA2~.C + oo)\. 

i>o |fc|<r2^/2] 

The sufficient conditions for the construction of shear let frames were given by [15j . 
Lemma 3.9. [Theorem 3.4 in USj/ J let (j),tp £ L^(M^) be functions such that 

0(6,6) < C2min{i,!eirnmin{i,ie2rn 

and 

1^(6 , 6) I < min{l, 16 r } min{l, 16 r^} min{l, 16 r^}, 

for some positive constants Ci,C2 < oo and a > 7 > 3. Define tp{xi,X2) = ^{x2,xi) and let 
Linf be defined by Li^f = ess infgg]g2 0(6 0). Suppose that there is a constant Li^f > such that 
< Linf ^ -^inf- Then there exist a sampling parameter c = (ci,C2) with ci = C2 stic/i that 
SH{(f), ijj, ijj] c) forms a frame for L'^(M?) 

Now we state the main results in this section. 

Theorem 3.10. For < a < 7r/2. Let aiq^{S,) , o^N^iCj ^he mask of the pseudo splines of type 
IL For I2 = 0, let Ni > N2 > 2. For I2 > 0, let Ni > N2 > 2. Let (^(6 be the associated refinable 
function with (iN^iO- L)efine the shearlet if) € L'^{E?) by 

i^{20 = Kii)m)m), e = (6,6)GK2, ^g^g^ 

where 6(6) = e~*^^a/vjj^(6 +71"). Then for given a € (0,7r/2), there exits a sampling constant 
c > such that the shearlet system '^{i[; c) forms a frame for {/ € L2(M^) : supp/ G Ci(a) UC3(a)} 
wit/i c = (ci, C2) G (K+)^ a'T'C^ C2 < ci < c. 

Proof. We first give the lower bound of the shearlet 'i/'(6- Let < a < 7r/2. By Theorem 13.41 we 
obtain 

|<A(6)I >C^X[-2«,2a](6) and |<^(6)l >C4X[-a,a](6)- (3-10) 

From (13. Sp . we have 

|6(6)l > |6(a)|x[- (6). (3.11) 

Combine (13.101) and (13.111) . we obtain 

m)\ = IK6)<^(6)<^(6)l > \b{a)\C'^Ci ■ xniO 
where = = (6)6), 6 ^ [—2a, —a] U [a, 2a], 6 ^ [—a, a]}. Since 

oo 

U U A2jSln = c, 

3=0 |fc|<r2^/2] 
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we conclude that 

= Yl mslA2-,0\' > UniXn{SlA2-,() > Uni on C, 

where U^i = {\b{a)\CiC4f . 

By p.Sp and ()3.8p . we obtain the upper bound of the shear let 

< Cmin{l, 1^1 • min{l, 1^11"'^^+"} • min{l, I^T'^^^+n- (3.12) 

Let /3jv^/ = 2N — k, for fixed N, /3jv,« decreases as / increases. For fixed /, fSjy^i increases as N 
increases. When I = N — 1, /^tv.z increases as N increase [9j. From the table 1, it easy to see that 

2Ni > 2N2 > 3, for l2 = 

and 

2Ni >2N2-K>3, for I2 > 0. 
Hence, the condition of Theorem 13.91 holds. Te theorem is proved. □ 

Remark 3.11. i&.y\) is a standard construction for separable wavelet basis of L2{M?) fBjj. Therefore, 
Theorem \3.10\ gives a connection between the shearlet analysis and wavelet analysis and also, we 
hope, enriches the theory of shearlet analysis. 

Since the decay rate of |i(?5>(OI is half of that of |20(OI- The results on pseudo splines of type I 
follow directly from Theorem 13.101 

Theorem 3.12. For < a < tt/2. Let cin^{C), o^n^{^) be the mask of the pseudo splines of type 
L For I2 = 0, let Ni > N2 > 3. For I2 > 0, let Ni > N2 > 8. Let 0(^) be the associated refinable 
function with ojq^j^i^)- Define a shearlet tp € L^(M^) by 

ij{2i) = S(ei)<^(6)<^(6), e = (6,6) G 

where 6(6) = e~*^^a7v^(6 + t^)- Then there exits a sampling constant c > such that the shearlet 
system "^(tp-jc) forms a frame for {/ € L2{M?) : supp/ G C{a)i Li €{0)3} for any sampling matrix 
Mc with c = (ci, C2) G (1^+)^ and C2 < ci < c. 

Remark 3.13. In 115], P. Kittipoom, G. Kutyniok, and W.-Q Lim gave a similar construction of 
the shearlet from pseudo splines of type I with order {N, I) such that I > 10 and ^ < N < 31 — 2. 
Compared with their results, our construction have smaller support. Moreover, we prove for any 
given < a < vr/2, one can construct shearlet frame. This gives more feasible choice in applications. 

4 Examples 

To illustrate our results, we give two examples in this section. 

Example 4.1. Let be the B- spline function of order 3. Define the shearlet 

i^i20 = 2-4e-*«Hl - e-'^^)''mi)m2) ^ = (6,6) G K^- 

Then for any given < a < 7r/2, there exits a sampling constant c > such that the shearlet 
system "^{ijj; c) forms a frame for {/ G L2{M?) : supp/ € Ci(a) U €3(0)} for C2 < ci < c. 
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Figure 1: i'{SkA22 -x), k = -2, -1,0,1,2, in Example O 

Example 4.2. Let B4 be the B-spline function of order 4. Define the shearlet 
M'^O = 2-^e-'^'il - e~'^'fM^i)B^i^2), ^ = (6,6) G M^. 

Let (j){x) = i?4(xi)i?4(x2) and '02(2;i,X2) = "01(^1) ^2)- Then there exit sampling constant c > 
such that the shearlet system ^'((/), -^i, -02! c) provides (almost) optimally sparse approximations of 
function f G <^'^{v) 
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